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THE  PATH  RESISTANCE  METHOD  FOR  BOUNDING  THE  SMALLEST  NONTRIVIAL 

EIGENVALUE  OF  A  LAPLACIAN* * * § 

STEPHEN  GUATTERYt,  TOM  LEIGHTON^,  AND  GARY  L.  MILLER§ 


Abstract.  We  introduce  the  path  resistance  method  for  lower  bounds  on  the  smallest  nontrivial  eigen¬ 
value  of  the  Laplacian  matrix  of  a  graph.  The  method  is  based  on  viewing  the  graph  in  terms  of  electrical 
circuits;  it  uses  clique  embeddings  to  produce  lower  bounds  on  A2  and  star  embeddings  to  produce  lower 
bounds  on  the  smallest  Rayleigh  quotient  when  there  is  a  zero  Dirichlet  boundary  condition.  The  method 
assigns  priorities  to  the  paths  in  the  embedding;  we  show  that,  for  an  unweighted  tree  T,  using  uniform 
priorities  for  a  clique  embedding  produces  a  lower  bound  on  A2  that  is  off  by  at  most  an  0(log  diameter(T)) 
factor.  We  show  that  the  best  bounds  this  method  can  produce  for  clique  embeddings  are  the  same  as  for  a 
related  method  that  uses  clique  embeddings  and  edge  lengths  to  produce  bounds. 

Key  words.  Laplacians,  Laplacian  spectra,  graph  eigenvalues  and  eigenvectors,  graph  embeddings 

Subject  classification.  Computer  Science,  Applied  and  Numerical  Mathematics 

1.  Introduction.  In  this  paper  we  consider  methods  based  on  graph  embeddings  for  estimating  the 
smallest  nontrivial  eigenvalue  of  the  Laplacian  matrix  representation  of  a  graph.  The  Laplacian  is  one  of 
many  ways  to  view  a  graph  as  a  matrix;  it  is  defined  as  follows:  Let  G  =  (V,  jF)  be  an  undirected  graph  with 
vertices  . . . ,  Vn-  Then  the  Laplacian  of  G  is  an  n  x  n  matrix  L  such  that 


hj  —  \ 


degree  (ui 


if  i  =  j 
if  (i,  j)  G  E 
otherwise 


It  is  not  hard  to  see  that  L  is  positive  semidefinite  (all  eigenvalues  are  >  0).  Since  the  row  sums  of  L  are  all 
zero,  the  smallest  eigenvalue  is  zero.  If  G  is  connected  the  second  smallest  eigenvalue  A2  is  positive. 

The  study  of  the  connection  between  Laplacian  spectra  (particularly  with  respect  to  A2)  and  properties 
of  the  associated  graphs  dates  back  to  Fiedler’s  work  in  the  1970’s  (see,  e.g.,  [Fie73]  and  [Fie75]).  These 
properties  have  been  used  in  graph  algorithms,  particularly  algorithms  for  finding  small  separators  [PSL90]. 

The  Laplacian  also  has  an  important  role  in  representing  physical  problems.  It  often  occurs  in  finite 
difference,  finite  element,  and  control  volume  representations  of  problems  involving  elliptic  partial  differential 
equations.  These  problems  often  include  a  Dirichlet  boundary  condition  that  specifies  that  the  values  at  a 
set  of  vertices  are  zero.  To  represent  this  condition  in  the  Laplacian,  the  rows  and  columns  corresponding  to 
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the  boundary  vertices  are  deleted  from  the  matrix.  The  resulting  matrix  is  positive  definite,  and  its  smallest 
eigenvalue  is  the  one  of  interest. 

Bounds  on  the  smallest  nonzero  eigenvalues  of  both  forms  of  Laplacian  have  other  important  applications. 
Since  the  matrices  are  .symmetric,  their  extreme  eigenvalues  can  be  used  in  computing  their  condition  num¬ 
bers.  which  are  used  in  the  study  of  iterative  linear  system  solvers  to  estimate  rates  of  convergence 
and  to  analyze  the  quality  of  preconditioners  [Axe92,  GMZ95].  Bounds  on  A2  are  useful  in  the  analysis 
of  spectral  partitioning,  both  because  A2  occurs  in  bounds  on  cut  quality  [Moh89],  and  because  they  can 
be  used  in  isolating  structural  properties  of  the  eigenvectors  used  in  making  the  cuts  [GM95.  ST96].  The 
eigenvalue  A2  has  been  related  to  expansion  properties  of  graphs,  and  can  be  used  in  determining  if  a  graph 
is  an  expander  [AM85a,  AI086]. 

Related  work  involving  graph  embeddings  has  been  used  to  bound  the  mixing  time  for  random  walks 
[JS89.  SJ89,  DS91,  Sin92,  Kah9r)].  In  this  case,  the  bound  is  on  the  second  largest  eigenvalue  of  the  Markov 
chain  transition  matrix.  Mixing  time  bounds  have  been  used  in  approximation  algorithms  for  a  number  of 
problems  (see,  e.g..  [DFK91]). 

In  many  of  the  preceding  applications,  it  is  necessary  to  show  bounds  for  classes  of  graphs  in  order  to 
state  results  in  terms  of  asymptotic  algorithm  behavior.  For  example,  in  the  analysis  of  spectral  partitioning, 
[GM95^  uses  eigenvalue  bounds  on  a  family  of  bounded-degree  graphs  to  prove  facts  about  the  structure  of 
eigenvectors  used  in  partitioning:  [ST96]  gives  an  upper  bound  on  planar  graph  eigenvalues  that  can  be 
applied  in  bounds  on  the  cut  quotient  of  the  resulting  cut.  The  embedding  techniques  we  present  are 
well-suited  to  producing  such  general  results,  and  can  be  used  with  known  results  about  embeddings.  We 
have  used  them  to  generate  lower  bounds  for  families  of  graphs  where  the  critical  path  resistance  can  be 
parameterized  in  terms  of  the  size  of  the  family  member. 

Note  that  the  applications  listed  above  involve  both  upper  and  lower  bounds.  Upper  bounds  on  A2  are 
usually  given  by  exhibiting  small  Rayleigh  quotients.  If  x  0  is  an  n-vector  its  Rayleigh  quotient  is  In 

fact  A2  =  min^^„^  Thus  a  R<iyleigh  quotient  for  a  vector  x  such  that  Y^:k~0  gives  an  upper  bound 

on  A2.  As  a  beautiful  example,  Spielman  and  Teng  have  proved  that  A2  =  0{l/n)  for  bounded  degree  planar 
graphs  by  showing  how  to  construct  a  vector  with  a  small  Rayleigh  quotient  for  each  such  graph  [ST96]. 

In  the  zero  Dirichlet  boundary  ca,se  we  are  interested  in  the  smallest  eigenvalue  A.  Consider  the  Laplacian 
L  with  the  rows  and  columns  of  the  boundary  vertices  deleted.  Then  A  =  minx  1^7^  •  In  this  case  the 
Rayleigh  quotient  has  a  continuous  analog:  Let  Q  be  a  compact  domain  in  IR"^  with  boundary  F,  e.g.  the 
unit  square  in  the  plane,  and  let  ?/  0  vary  over  all  differentiable  functions  defined  on  Q  such  that  ii  =  0  on 

r.  Then  the  Rayleigh  quotient  is 


c  = 


min 

ii 


where  Vu  is  the  gradient  of  n.  This  lower  bound  on  c  is  known  as  a  Poincare  inequality.  This  analog  provides 
useful  insight  on  the  problem  and  helps  with  the  proofs  to  follow. 

We  start  by  describing  a  very  simple  yet  powerful  method  for  bounding  A2  from  below  (Figure  1).  We 
show  that  the  Path  Resistance  Method  correctly  bounds  A2  from  below.  Note  that  the  quality  of  the  lower 
bound  improves  with  the  cpiality  of  the  path  embedding  since  better  embeddings  decrease  congestion.  We 
also  show  that  it  is  sometimes  possible  to  improve  the  lower  bound  by  assigning  a  priority  to  each  path  and 
allocating  the  edge  conductances  ba.sed  on  the  priorities.  The  path  resistance  method  as  stated  in  the  figure 
does  not  mention  priorities:  we  refer  to  that  version  as  the  unit-priority  case. 
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Procedure:  The  Path  Resistance  Method 
Method: 

1.  Construct  a  path  embedding  of  Kn  into  i.e., 
construct  a  path  between  every  pair  of  vertices  of 
G. 

2.  For  each  edge  eij  in  G  compute  the  congestion 

^ij  —  ^  ^}l 

3.  For  each  path  P  and  each  edge  eij  on  P  allocate 
a  resistor  of  size  Cij  to  P. 

4.  For  each  path  P  compute  its  resistance,  i.e., 

X)e  eP  ^  maximum  resistance  over 

all  paths. 

5.  return  “n/r  <  A2”. 

Fig.  1.1.  The  Uniform  Priority  Path  Resistance  Method 

The  clique  embedding  version  of  the  Path  Resistance  Method  is  closely  related  to  another  clique  embed¬ 
ding  technique  presented  in  [Kah96]  that  assigns  edge  lengths  and  computes  the  sum  of  the  lengths  of  all 
paths  incident  to  each  edge.  We  show  that  these  methods  are  duals  of  one  another,  and  that  the  best  lower 
bounds  computed  by  these  methods  are  the  same. 

We  also  present  a  version  of  the  Path  Resistance  Method  that  is  applicable  in  the  zero  Dirichlet  boundary 
case.  Instead  of  embedding  a  clique,  we  embed  a  star  into  the  graph,  with  the  boundary  vertices  consolidated 
into  the  vertex  at  the  center  of  the  star. 

Finally,  we  show  that  the  chque  embedding  version  of  the  Path  Resistance  Method  with  uniform  priorities 
is  robust  when  applied  to  trees.  We  produce  upper  and  lower  bounds  to  show  that  for  a  tree  T  using  uniform 
priorities,  our  estimate  is  off  by  at  most  a  factor  of  O (log  diameter (T)).  It  is  somewhat  surprising  that  such 
a  simple  way  of  setting  the  priorities  gives  such  a  good  bound.  It  is  open  as  to  how  well  the  uniform  priority 
method  does  for  general  graphs. 

2.  Previous  Work.  Prior  work  related  to  the  applications  of  the  lower  bounds  has  been  cited  in  the 
Introduction.  The  use  of  clique  embeddings  to  bound  eigenvalues  arose  in  the  analysis  of  mixing  times  for 
Markov  chains  by  Jerrum  and  Sinclair  [JS89]  [SJ89].  Further  work  in  this  direction  was  done  by  Diaconis 
and  Strook  [DS91]  and  by  Sinclair  [Sin92].  Kahale  [Kah96]  generalized  this  work  in  terms  of  methods 
assigning  lengths  to  the  graph  edges,  and  showed  that  the  best  bound  over  all  edge  length  assignments  is 
the  largest  eigenvalue  of  the  matrix  P^F,  where  F  is  a  matrix  representing  the  path  embedding  ([Kah96] 
also  cites  unpublished  work  by  Fill  and  Sokal  in  these  directions).  He  also  gave  a  semidefinite  programming 
formulation  for  a  model  allowing  fractional  paths,  and  showed  that  the  bound  is  off  by  at  most  a  factor  of 
log^  n.  He  proved  this  gap  is  tight;  he  also  noted  that  the  results  can  be  applied  to  Laplacians  with  suitable 
modifications. 

3.  Terminology,  Notation,  and  Background  Results.  We  assume  that  the  reader  is  familiar  with 
the  basic  definitions  of  graph  theory  (in  particular,  for  undirected  graphs),  and  with  the  basic  definitions 
and  results  of  matrix  theory.  A  graph  consists  of  a  set  of  vertices  V  and  a  set  of  edges  E\  we  denote  the 
vertices  (respectively  edges)  of  a  particular  graph  G  as  V(G)  (respectively  E{G))  if  there  is  any  ambiguity 
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about  which  graph  is  referred  to.  When  it  is  clear  which  graph  we  arc  referring  to,  we  use  r?  to  denote  \V\. 

We  use  the  term  path  graph  for  a  tree  that  has  exactly  two  vertices  of  degree  one.  That  is,  a  path 
graph  is  a  graph  consisting  of  exactly  its  maximal  path.  A  star  is  a  tree  with  exactly  one  vertex  that  is  not 
a  leaf.  We  call  the  non-leaf  vertex  the  center  of  the  star. 


3.1.  Matrices  and  Matrix  Notation.  Wc  use  capital  letters  to  represent  matrices  and  bold  lower¬ 
case  letters  for  vectors.  For  a  matrix  A,  or  [A],^  represents  the  clement  in  row  i  and  column  j:  for  the 
vector  X,  ,r,  or  [xj,  represents  the  entry  in  the  vector.  The  notation  x  =  0  indicates  that  all  entries  of  the 
vector  X  are  zero;  1  indicates  the  vector  that  has  1  for  every  entry.  For  Laplacians,  we  index  the  eigenvalues 
of  an  77  X  77  matrix  in  non-decreasing  order:  Aj  represents  the  smallest  eigenvalue,  and  A^  the  largest.  We 
use  the  notation  Ai(L)  (respectively  Xi{G))  to  indicate  the  eigenvalue  of  matrix  L  (respectively  of  the 
Laplacian  of  graph  G)  if  there  is  any  ambiguity  about  which  matrix  (respectively  graph)  the  eigenvalue 
belongs  to.  u,  repre.sents  tlie  eigenvector  corresponding  to  A,-. 


3.2.  The  Laplacian  Matrix  Representation  of  a  Graph.  A  common  matrix  representation  of 
graphs  is  the  Laplacian.  Let  D  be  the  matrix  with  =  dcgrce(7;,)  for  Vi  G  V{G),  and  all  off-diagonal 
entries  equal  to  zero.  Let  A  be  the  adjacency  matrix  for  G  [a-ij  =  1  iff  G  E{G).  0  otherwise).  Then 

the  Laplacian  representation  of  G  is  the  matrix  L  —  D  ~  A. 

The  following  are  some  useful  facts  about  the  Laplacian  matrix: 

•  The  Laplacian  is  symmetric  positive  semidefinite,  so  all  its  eigenvalues  are  greater  than  or  equal  to 
0  (see  e.g.  [AM85b]). 

•  A  graph  G  is  connected  if  and  only  if  0  is  a  simple  eigenvalue  of  its  Laplacian  (see  e.g.  [AM85b]). 

•  The  following  characterization  of  A2  holds  (see  e.g.  [Fie73]): 


(.3.1) 


.  x'^  Lx 

A2  ==  mm  — ~ — . 

xJ.r  X^x 


•  For  any  vector  x  and  Laplacian  L  of  the  graph  G,  we  have  (see  e.g.  [Fie75]): 


(3.2)  Lx  =  ^  (.r,  -  Xj)^ 

(r..r,)eE((7) 

An  edge- weighted  graph  has  a  real,  nonzero  weight  vuj  associated  with  each  edge  (i’,.  t’j)  (we  consider 
a  zero  edge  weight  to  indicate  the  lack  of  an  edge).  Fiedler  extended  the  notion  of  the  Laplacian  to  graphs 
with  positive  edge  weights  [Fie75];  he  referred  to  this  representation  as  the  generalized  Laplacian.  Let 
be  the  (positive)  weight  of  edge  (ij)  in  graph  G.  Then  the  entries  of  the  generalized  Laplacian  L  of 
G  are  defined  as  follows:  is  the  sum  of  the  weights  of  the  edges  incident  to  vertex  Vi:  for  ?'  ^  j  and 

(r, ,  Vj)  G  E{G).  Uj  —  —v'jj,  and  /,j  —  0  otherwise. 

With  the  exception  of  equation  (3.2),  the  properties  listed  above  also  apply  to  generalized  Laplacians. 
A  slightly  modified  version  of  (3.2)  holds  for  the  generalized  Laplacian  L: 


x^Lx=  Y, 

(v..v,)CE(0) 

We  usually  use  L  to  denote  the  Laplacian  of  a  graph  G,  and  K  to  represent  the  Laplacian  of  the  clique 
K  (whether  K  refers  to  the  graph  or  the  Laplacian  will  be  clear  from  context). 

3.3.  Graph  Embeddings.  We  start  with  the  notion  of  embedding  one  graph  into  another.  Let  G  and 
H  be  connected  graphs  such  that  the  vertex  set  of  H  is  a  siibset  of  the  vertex  set  of  G.  An  embedding 
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of  H  into  G  is  a  collection  F  of  path  subgraphs  of  G  such  that  for  each  edge  {vi^Vj)  €  E{H)^  F  contains  a 
simple  path  7^j  from  Vi  to  Vj  in  G.  For  full  generality,  we  allow  fractional  paths  in  our  embeddings:  i.e., 
an  edge  {vi^Vj)  €  E{H)  can  be  associated  with  a  finite  collection  of  simple  paths  from  Vi  to  Vj  in  G\  each 
such  path  has  a  positive  fractional  weight  associated  with  it  such  that  the  weights  add  up  to  1.  If  a  path  7 
includes  edge  e,  we  say  that  7  is  incident  to  e. 

In  the  unit-priority  case,  the  congestion  cong(e)  of  an  edge  e  G  E(G)  is  the  number  of  paths  in  the 
embedding  incident  to  e  (fractional  paths  using  e  contribute  their  fractional  weight  to  this  count).  The 
resistance  r(e)  of  an  edge  is  its  congestion  divided  by  its  weight:  r(e)  =  The  resistance  of  a  path 

r(7)  is  computed  from  the  edge  resistances  in  the  same  way  that  the  resistance  of  a  series  of  resistors  would 
be  computed  in  an  electrical  network.  That  is,  r(7)  =  ^(e).  The  largest  path  resistance  r^ax  is  defined 

as  Tmax  =  max^^p  ^(7)-  The  corresponding  path  is  called  the  critical  path.  Note  that  r^ax  is  always 
defined  with  respect  to  an  embedding;  the  particular  embedding  will  be  clear  from  the  context,  and  thus  we 
will  not  specify  it  in  the  notation. 

The  definitions  change  slightly  when  we  use  path  priorities.  Let  pij  be  the  path  priority  for  7^^,  and  let  p 
be  the  vector  of  priorities.  The  congestion  cong(e)  of  edge  e  is  the  sum  of  the  priorities  of  e’s  incident  paths. 
Each  edge  has  a  resistance  factor  r(e),  defined  as  the  congestion  of  the  edge  divided  by  its  weight.  The 
resistance  of  edge  e  with  respect  to  ^ij  is  The  resistance  r(7)  of  a  path  7  is  the  sum  of  the  resistance 
factors  of  the  edges  in  the  path  divided  by  the  path  priority  (i.e.,  r{'yij)  =  ^  ^(^))*  The  maximum 

path  resistance  taken  over  all  paths  in  F  is  denoted  fmax-  The  corresponding  path  is  called  the  critical 
path. 

We  note  that  using  these  modified  definitions,  the  general  method  for  computing  lower  bounds  using 
priorities  is  clear:  Compute  rmax  according  to  the  revised  definitions.  Then  <  A2. 

4.  Lower  Bounds  on  A2  for  Generalized  Laplacians. 

4.1.  A  Bound  Based  on  a  Clique  Embedding.  We  now  consider  the  case  in  which  G  is  a  graph  for 
which  we  would  like  to  bound 'A2  of  G’s  (generalized)  Laplacian  from  below.  We  will  use  a  clique  embedding 
to  decompose  the  graph.  The  decomposed  graph  is  then  viewed  in  terms  of  an  electrical  circuit  analogy. 

We  start  with  an  easy  lemma  that  helps  convey  the  nature  of  the  bounding  technique.  Let  L  be  the 
Laplacian  of  (edge- weighted)  graph  G,  and  let  K  be  the  Laplacian  of  the  complete  graph  on  n  =  |V(G)| 
vertices. 

Lemma  4.1.  Let  a  be  a  positive  real  number.  If  aL  —  K  is  positive  semidefinite,  then  \2{L)  > 

Proof  Assume  that  aL  —  K  is  positive  semidefinite.  Consider  U2  (aL  — iF)u2,  where  U2  is  the  eigenvector 
corresponding  to  A2(L).  Without  loss  of  generality,  assume  that  U2  has  been  normalized  to  length  1.  Note 
that  K  :=  nl  ~  J,  where  J  is  the  n  x  n  matrix  with  all  entries  equal  to  1.  Thus  it  is  clear  that  every  vector 
orthogonal  to  1  is  an  eigenvector  of  K  with  eigenvalue  n.  Therefore 

U2  (aL  —  K)u2  =  auj Lu2  -  U2  iFu2  =  aA2  —  n  >  0, 

where  the  last  inequality  holds  by  the  assumption  that  aK  —  L  is  positive  semidefinite.  The  lemma  follows. 

□ 

Thus,  if  we  can  demonstrate  an  a  sufficiently  large  to  make  aL  —  K  positive  semidefinite,  we  can  bound 
A2  from  below.  To  do  this,  we  form  an  embedding  (F,  p)  of  the  complete  graph  into  G,  then  break  aL—K  into 
pieces  consisting  of  edges  in  K  and  their  corresponding  paths  in  F.  We  use  the  following  additional  notation: 
Eij  denotes  the  Laplacian  of  the  graph  on  V{G)  that  has  only  the  single  edge  (i,i).  L{'yij)  represents  the 
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gcnrralizrd  Laplacian  on  V((7)  whose  edges  are  the  edges  of  and  whose  edge  weights  are  set  as  follows: 
each  path  "'/,j  using  an  edge  (<7,  h)  gets  a  share  of  proportional  to  ■  This  works  out  to 

Ptj  ^  Pu 
cong{gJi)  r^k  ' 

It  is  easy  to  see  that  the  weights  assigned  to  all  incident  paths  for  {g.  h)  sum  to  Wgh,  (If  fractional  paths  arc 
involved,  then  the  matrix  L(7,j)  is  multiplied  by  the  fractional  coefficient  of  path 
\\c  now  state  the  theorem  that  gives  us  the  technique  for  bounding  A2: 

Thf:orem  4.2.  For  any  generalized  LaplaeAan  L  of  a  connected  graph  G  xiuth  positive  edge  weights,  and 
any  ehqnr  embedding  into  G, 


Enmx 


Proof 

Any  embedding  of  the  complete  graph  into  G  defines  a  natural  decomposition  of  the  Laplacian  L  in 
terms  of  the  fj's:  define  the  matrix  Lj^nnpnfh  as  the  Laplacian  of  the  graph  on  V{G)  that  includes  any  edges 
of  G  not  used  in  any  path  It  is  easy  to  see  that 

L  ^  ^  T(q)  4*  Fnonpnt h  ‘ 

7cr 

Since  Lj,onpnth  is  positive  semidefinite,  the  following  inequality  holds  for  any  positive  a  and  real  vector  x: 


(4.1) 


X  <  (aL  -  A")  X. 


i<J 


Applying  properties  of  linearity  and  rewriting  gives 


(•L2)  (ox^L(q,j)x  -  x^Ejjx)  <  x^  (aL  ~  A")  x. 

Now  consider  the  terms  ox^  L(q,j)x  -  x^Eijx.  By  (3.2),  x^E.jx  =  (.?;,  -  Likewise,  the  term  for  the 

path  can  be  written  in  terms  of  its  edges  and  edge  weights: 


ax^  L{')ij)x  =  a  Y]  — (x, 

Since  all  the  priorities  and  resistance  factors  are  positive,  the  sum  in  the  expression  above  can  be  interpreted 
as  the  energy  dissipation  of  a  series  of  resistances  (the  reciprocals  of  the  ^’s)  given  a  set  of  potentials  (the 
T,  s)  at  the  nodes  between  the  resistances.  It  is  well  known  (see  e.g.  [DS84])  that  this  quantity  is  minimized 
when  the  potentials  at  the  internal  nodes  are  consistent  with  Kirchoff’s  law  (this  is  easily  seen  through  an 
application  of  the  Cauchy-Schwarz  inequality).  At  that  minimum  the  sum  is  the  square  of  the  potential 
differences  of  the  endpoints  divided  by  the  path  re,sistance: 


Thus,  if  a  =  we  have  for  every  path  7,^  and  every  x  that 

ax^  L(7,j)x  -  x^E.jX  >  -  (.t,  -  xjf  >  0. 
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For  this  value  of  a,  the  left-hand  sides  of  inequalities  (4.1)  and  (4.2)  are  nonnegative.  Thus  r^ax^  -  K  is 
positive  semidefinite,  and  the  theorem  holds  by  Lemma  4.1. 

If  fractional  paths  are  involved,  then  each  L{'yij)  is  multiplied  by  the  fractional  coefficient,  as  is  the 
corresponding  Eij.  The  fractional  coefficients  cancel,  so  it  is  easy  to  see  that  the  theorem  also  holds  in  this 
case.  □ 

We  now  compare  the  path  resistance  method  with  the  edge  length  method  described  in  [Kah96].  The 
method  as  applied  to  Laplacians  is  described  as  follows: 

•  Specify  a  clique  embedding  for  the  graph,  and  assign  each  edge  a  positive  length. 

•  Compute  the  length  of  each  path  with  respect  to  the  edge  lengths. 

•  For  each  edge,  compute  the  sum  of  the  lengths  of  all  incident  paths  divided  by  the  length  of  that 
edge.  Let  Pm^x  be  the  maximum  such  value  taken  over  all  the  edges.  Then  is  a  lower  bound 
for  A2. 

The  proof  that  this  produces  a  lower  bound  can  be  done  in  a  fashion  similar  to  the  proof  given  above 
for  the  path  resistance  method;  however,  the  argument  manipulates  the  x^i^x  term.  More  specifically,  the 
difference  across  each  edge  (ij)  €  K  gets  written  as  a  telescoping  sum  of  the  differences  along  7^^.  The 
Cauchy-Schwarz  inequality  is  applied,  and  the  resulting  sum  is  reorganized  in  terms  of  edges  in  L. 

In  the  unit  priority/unit  length  case,  the  two  methods  can  be  thought  of  in  terms  of  the  following  electrical 
analogies:  the  path  resistance  method  partitions  the  conductance  (i.e.,  the  reciprocal  of  the  resistance)  of 
an  edges  e  G  E{G)  equally  among  its  incident  paths.  The  resulting  paths  have  path  conductances;  we  find 
a  multiplier  (rmax)  sufficient  to  increase  the  conductance  of  every  path  so  that  it  is  at  least  as  big  as  the 
conductance  of  the  corresponding  clique  edge.  The  edge  length  method  also  deals  with  path  conductances, 
but  instead  assigns  each  path  a  demand  equal  to  the  path’s  length;  this  insures  that  each  path  has  unit 
conductance  to  support  a  clique  edge.  We  find  a  multiplier  (/9max)  sufficient  to  make  each  edge’s  conductance 
at  least  as  large  as  the  sum  of  the  demands  of  its  incident  paths.  In  both  methods,  the  multiplier  shows  up 
in  the  bound. 

The  two  methods  are  duals  of  each  other  in  the  following  sense: 

Theorem  4.3.  Given  a  path  embedding,  the  best  lower  bound  for  the  path  resistance  method  taken  over 
all  allowed  priority  assignments  and  the  best  lower  bound  for  the  edge  length  method  taken  over  all  allowed 
length  assignments  are  the  same. 

Proof  To  show  this,  we  will  use  a  representation  presented  by  Kahale  in  [Kah96].  He  defines  an 
embedding  matrix  F  as  follows:  Each  row  of  F  represents  a  path  between  a  pair  of  distinct  vertices.  A  row 
entry  is  1  if  the  corresponding  edge  is  in  the  path  and  zero  otherwise  (this  is  easily  generalized  to  fractional 
paths).  It  is  easy  to  see  that  multiplying  a  vector  of  edge  lengths  by  F  computes  the  length  of  each  path,  and 
that  multiplying  the  result  by  F^  sums  the  lengths  of  the  incident  paths  for  each  edge.  Thus  multiplying  an 
edge  vector  by  F^F  and  dividing  the  result  for  each  edge  by  the  original  length  and  taking  the  maximum  is 
a  way  to  compute  the  value  pmax  Rsed  in  the  edge  length  lower  bound.  F^F  is  clearly  a  nonnegative  matrix, 
and  Kahale  shows,  using  the  properties  of  such  matrices  (see,  e.g.,  [Min88]),  that  the  best  lower  bound  for 
this  method  is  the  largest  eigenvalue  of  this  matrix. 

The  path  resistance  lower  bound  can  be  computed  by  multiplying  a  vector  of  priorities  by  FF^,  then 
dividing  the  result  termwise  by  the  original  priorities  and  setting  rmax  to  the  maximum  resulting  value.  The 
same  arguments  about  nonnegative  matrices  apply.  However,  FF^  and  F^F  have  the  same  set  of  nonzero 
eigenvalues;  this  is  easily  seen  by  showing  that  for  every  eigenvector  u  of  a  nonzero  eigenvalue  of  FF^,  F^u 
is  an  eigenvector  of  F^F  with  nonzero  eigenvalue;  likewise,  if  u  is  an  eigenvector  of  a  nonzero  eigenvalue  of 
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r^r.  Fu  is  an  eigenvector  of  FF^  with  nonzero  eigenvalue.  Thus  the  best  lower  bounds  produced  arc  the 
same.  □ 

4.2.  Star  Embeddings.  Embeddings  other  than  clique  embeddings  are  useful.  In  particular,  we  embed 
stars  into  graphs  with  zero  Dirichlct  boundary  conditions.  Such  boundary  conditions  specify  a  set  of  vertices 
whose  values  are  zero.  Let  G  be  a  graph  with  such  a  zero  boundary,  and  let  L  be  its  Laplacian.  Let  X  be 
the  set  of  non-zero  vectors  consistent  with  the  boundary  restriction.  We  are  interested  in  that  value  of  the 
following  quantity: 

.  x^Lx 

A  =  mm  — vp — 

xC  A'  X 

To  find  a  lower  bound,  we  embed  a  star  into  G  such  that  every  boundary  vertex  is  mapped  to  the  center 
of  the  star,  and  every  other  vertex  is  mapped  to  a  distinct  leaf.  Thus  we  specify  a  path  in  G  from  each 
non-boundary  vertex  to  a  boundary  vertex.  Priorities,  congestions,  resistances,  and  r^ax  are  all  defined  as 
before.  The  following  theorem  applies: 

Theorem  4.4. 


For  brevity,  we  omit  the  proof,  which  is  almost  the  same  as  the  proof  in  the  clique  case.  For  readers  wishing 
to  construct  the  proof  themselves,  we  note  that  for  the  star,  the  Rayleigh  quotient  of  any  nonzero  vector 
with  zero  at  the  center  is  exactly  1. 

Matrices  for  the  zero  Dirichlct  boundary  case  often  represent  boundary  vertices  implicitly.  That  is.  the 
rows  and  columns  of  boundary  vertices  are  deleted  from  the  Laplacian  of  the  full  graph  G.  In  that  cavsc, 
edges  between  non-boundary  and  boundary  vertices  are  not  explicitly  represented  by  off-diagonal  entries  in 
the  resulting  matrix.  Instead,  they  show  up  as  surpluses  in  the  diagonal  entries  for  the  non-boundary  end 
of  the  edge. 

4.3.  Paths  with  Masses.  To  simplify  the  computation  of  quantities  such  as  congestions  for  embed¬ 
dings  into  trees  in  the  unit  priority  case,  we  introduce  the  idea  of  mass.  Consider  a  tree  T  and  a  unit  priority 
embedding  of  a  graph  H  into  into  T  (note  that  this  embedding  F  is  unique).  Using  the  definitions  from 
Section  3.3  above,  we  can  compute  the  critical  path.  We  distinguish  two  types  of  vertices  that  do  not  lie  on 
the  critical  path:  those  not  in  any  path  containing  a  critical  path  edge  (noncontributing  vertices),  and  those 
that  are  in  some  such  path  (contributing  vertices).  For  each  contributing  vertex  u,  there  is  a  first  critical 
path  vertex  tv  on  any  path  that  includes  both  v  and  a  critical  path  edge.  Note  that  since  we  are  working 
with  a  tree.  t\.  is  the  same  for  every  such  path.  We  will  call  tv  the  connection  point  for  v.  Assume  that 
the  critical  path  has  k  vertices  indexed  from  1  to  A*  in  order.  We  assign  ma.s.ses  to  the  critical  path  vertices 
as  follows:  Vertex  u,  has  a  mass  equal  to  1  plus  the  number  of  contributing  vertices  that  have  as  a 
connection  point. 

We  refer  to  the  graph  P  consisting  of  a  path  on  k  vertices  plus  the  masses  for  the  critical  path  as  the 
path-with-masses  model.  It  has  an  embedding  of  a  subgraph  of  H  consisting  of  a  subset  of  F  that  includes 
paths  with  both  endpoints  on  the  critical  path.  Congestions  are  computed  as  before,  except  that  now  each 
path  contributes  an  amount  equal  to  the  product  of  the  masses  of  its  endpoints  to  the  edge  congestions.  It 
is  easy  to  see  that  the  path  congestion  of  P  in  this  model  is  equal  to  the  conge.stion  on  the  critical  path  of 
the  original  embedding. 
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We  can  also  use  the  path-with-masses  model  in  computing  upper  bounds.  Let  y  be  a  vector  of  length 
k.  Let  L{P)  be  the  Laplacian  of  P,  and  let  M  be  the  diagonal  matrix  with  entry  \M]u  =  mi-  Let  L  be 
the  Laplacian  of  T,  the  tree  we  are  working  with.  Construct  a  vector  x  on  n  vertices  as  follows:  For  the  z* 
vertex  on  the  critical  path,  set  the  corresponding  entry  of  x  to  y,.  For  each  noncontributing  vertex,  set  the 
corresponding  entry  of  x  to  0.  For  each  contributing  vertex,  set  the  corresponding  entry  of  x  to  the  value  of 
y  for  its  connection  point.  Recall  Laplacian  Property  (3.2)  from  Section  3.2,  and  note  that  the  only  edges 
of  T  with  nonzero  differences  across  them  lie  on  the  critical  path.  Note  also  that  the  number  of  vertices  in 
T  assigned  the  value  yi  is  equal  to  rrii.  It  is  therefore  easy  to  see  that 

x^Lx  _  y^L(P)y 
x^x  y^My 

In  the  star  case  with  Dirichlet  boundary,  it  is  easy  to  see  that  if  we  enforce  the  boundary  condition  in  y, 
then  it  is  enforced  in  x,  which  imphes  that  is  an  upper  bound  on  A  as  defined  above.  For  clique 

embeddings,  note  that  enforcing  the  condition  y^Ml  =  0  is  equivalent  to  enforcing  the  condition  x^l  =  0; 
in  this  case,  ^yT^y^  is  an  upper  bound  on  A2. 

5.  Lower  Bounds  for  Trees.  In  this  section  we  show  that,  for  an  arbitrary  unweighted  tree  T  with 
diameter  diam(r),  assigning  all  paths  priority  1  in  the  clique  embedding  gives  a  lower  bound  within  a 
logdiam(r)  factor  of  A2.  The  argument  works  by  first  showing  the  bound  holds  for  a  tree  with  a  single  zero 
boundary  vertex,  then  applying  this  fact  along  the  critical  path  to  handle  the  chque-embedding  case. 

We  start  with  a  Dirichlet  boundary  case,  where  A  is  defined  as  in  Section  4.2.  Let  k  be  the  length  of  the 
critical  path. 

Theorem  5.1.  For  any  unweighted  tree  with  a  single  zero  boundary  vertex,  let  r^ax  be  the  maximum 
path  resistance  for  the  star  embedding  in  the  unit  priority  case.  Then 

max(lo^^^(^) 

^max 


Proof.  Note  that  if  the  boundary  vertex  separates  the  tree  into  multiple  components,  the  components 
can  be  considered  separately.  We  need  only  consider  the  component  with  the  critical  path. 

Let  L  be  the  Laplacian  of  the  tree.  For  any  such  tree  and  boundary,  we  construct  a  nonzero  vector  x 
consistent  with  the  boundary  condition  such  that 


(5.1) 


x^I/X  ^  c  max  (log  k,  1) 


where  c  is  a  constant  independent  of  the  choice  of  tree  and  boundary.  By  the  definition  of  A,  the  theorem 
then  holds.  The  construction  of  x  is  done  with  respect  to  the  path-with-masses  model  of  the  critical  path 
as  described  in  Section  4.3. 

Let  P  be  the  critical  path.  By  assumption  it  has  k  edges.  We  index  its  vertices  as  follows:  vertex  0 
is  the  boundary  vertex.  Vertex  i  is  the  vertex  at  distance  i  from  the  boundary.  The  maximum  index  is  k. 
Mass  rui  is  assigned  to  Vi  as  specified  in  Section  4.3  to  account  for  paths  that  start  off  the  critical  path,  but 
intersect  with  it. 

Because  we  are  dealing  with  the  critical  path  plus  a  zero  boundary,  we  can  express  the  value  Xi  as 
follows:  Let  6i  =  Xi  ~  Xi-i.  Then  Xi  =  allows  us  to  rewrite  some  of  the  quantities  from 

Inequality  (5.1)  as  follows: 

1  _  1  _  1 

(^*^)  •  ’ 

rmax  Z2j^i  ^3  zZi^l  ^ 
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and 


x'Lx 


E* 


Wc  specify  X  by  sottiiif^^  S,  =  i~~  2 .  For  A*  =  1  it  is  clear  that  the  quantities  on  the  right-hand  sides  of  (5.2) 
and  (5.3)  are  equal  and  the  theorem  holds.  For  A'  >  1,  we  have 

t-l  jr-l 

Tims  if  we  show  that  the  quantities  in  the  denominators  of  the  right-hand  sides  of  (5.2)  and  (5.3)  arc  within 
a  constant  factor,  the  theorem  holds.  Note  that  the  two  denominators  can  be  thought  of  as  expressions  in 
the  mys.  so  it  will  suffice  to  show  that  the  corresponding  coefficients  for  each  m,  are  within  a  constant  of 
each  other.  Thus,  we  want  to  show  that  i  is  within  a  constant  factor  of  .  For  i  =  1  these  arc 

both  1.  For  ?  >  1,  the  summation  can  be  bounded  using  integration  techniques;  in  this  case  the  sum  lies 
between  ?2  and  2/2.  Thus  the  square  of  the  sum  lies  between  i  and  4?,  which  gives  the  desired  result.  □ 

It  is  po.ssible  to  find  examples  for  which  this  gap  is  tight.  Consider  a  tree  consisting  of  a  path  graph 
connecting  the  center  vertices  of  a  series  of  stars.  In  particular,  let  the  path  have  A'  -f  1  vertices  numbered 
from  0  to  A'.  Vertex  r,  (0  <  ?  <  A*)  on  the  path  is  the  center  of  a  star  with  -  1  leaves.  For  simplicity 

of  reference,  we  refer  to  this  tree  as  the  bad  tree  for  parameter  k.  Vertex  uq  is  the  zero  Dirichlct  boundary. 

We  show  that  the  gap  for  the  bad  tree  is  tight  by  showing  a  set  of  priorities  that  give  a  lower  bound 
within  a  constant  of  the  upper  bound  constructed  as  per  the  proof  of  Theorem  4.4.  Let  T  be  a  bad  tree  with 
A'  >  1:  let  L  be  the  Laplacian  of  T.  The  path  connecting  the  stars  is  obviously  the  critical  path.  W’e  assume 
that  the  path  vertices  are  numbered  as  in  the  definition  of  tlie  bad  tree.  Path  edges  arc  numbered  by  the 
larger  of  the  endpoint  indices:  i.e.,  edge  i  is  from  Vj_.i  to  Vj. 

Theoricm  5.2.  Lrf  ro  6c  the  zero  houndari/  r)erf.cx  for  bad  tree  T  with  parameter  A*  >  1.  Let  f/jg 

unit  priority  lower  bound  estimate  on  X  of  L  given  that  boundary.  Then  A  =  0 

Proof  W'e  again  use  the  “path  witli  masses”  model.  It  is  clear  that  vertex  Vi  has  mass  We  note 

that  the  quantity  ^  is  always  greater  than  or  equal  to  1  in  the  allowed  range,  so  its  ceiling  is  bigger  by  no 
more  than  a  factor  of  2.  Because  we  are  interested  in  bounds,  we  drop  the  ceilings;  it  is  easy  to  verify  that 
this  does  not  change  the  results  below  by  more  than  a  constant  factor. 

As  in  the  proof  of  Theorem  4.4.  we  can  construct  a  vector  x  with  ,to  =  0  and,  for  z  >  0,  .  Values 

at  vertices  off  the  critical  path  are  set  as  per  Section  4.3,  Consider  the  Rayleigh  quotient  As  in  the 

previous  proof,  the  numerator  is  0(logA').  The  denominator  is 


=E-- 


As  noted  in  the  proof  to  Theorem  4.4,  3  ^  is  v/?  to  within  a  constant  factor,  so,  to  within  a  constant 

factor  the  denominator  is 

^  ^  7,2  ^ 

7T  =  E  T  ^  e(A‘"logA’) 


(note  that  k  >  1  by  assumption).  Thus  the  upper  bound  is  B(p). 
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For  the  unit  priority  lower  bound,  it  is  obvious  that  the  critical  path  runs  from  the  boundary  to  vertex 
Ufc.  As  noted  previously,  we  can  write  the  formula  for  r^ax  as  follows  (once  again  we  drop  the  ceilings  without 
changing  the  result  by  more  than  a  constant  factor): 


=  =  =e(fe^logfe). 

i=l  i=l 

Thus  the  upper  and  lower  bounds  differ  by  a  logfc  factor,  the  value  of  the  Rayleigh  quotient  numerator. 

We  can  get  a  better  lower  bound  and  close  the  gap  by  assigning  priorities  to  paths.  The  relative  sizes  of 
the  priorities  determine  the  relative  shares  of  the  conductance  that  the  paths  get.  In  the  current  problem, 
the  resistance  of  the  critical  path  is  high  because  it  gets  very  small  shares  of  the  conductances  of  the  edges 
close  to  the  boundary.  By  increasing  the  shares  of  longer  paths,  we  can  offset  this  problem;  because  the 
number  of  shorter  paths  using  edge  i  grows  fast  as  the  boundary  is  approached,  the  number  of  paths  getting 
an  increased  share  of  conductance  is  relatively  small  and  doesn’t  increase  the  resistance  of  these  shorter 
paths  too  much. 

We  now  specify  a  set  of  priorities  that  gives  us  a  lower  bound  that  matches  the  upper  bound.  Each  path 
that  either  starts  at  Vi  on  the  critical  path,  or  that  has  v,  as  its  connection  point,  gets  priority  Pi  =  VI. 
With  this  change,  the  congestion  of  edge  i  is  equal  to 

k  ^  7,2  ^ 

j=i  j=i  3=^ 

For  all  i  greater  than  1,  we  use  integration  techniques  to  get  the  bounds 


(5.4) 


2  >  2 

j=i 


1 


and 

(5.5) 


2  <  2/c^ 


For  z  =  1  the  upper  bound  is  replaced  by  2  ^2  -  the  lower  bound  remains  unchanged. 

Recall  that  in  the  case  where  path  7  has  priority  p,  we  compute  the  path  resistance  as  follows: 

^(7)  =  “y^cong(e). 

^  eG7 

Plugging  in  the  upper  and  lower  bounds  on  edge  congestion  in  (5.4)  and  (5.5)  above  gives  upper  and  lower 
bounds  on  the  path  resistance  for  a  path  7^  starting  at  Vj,  We  start  with  the  lower  bound;  the  third  line 
below  follows  from  an  application  of  standard  integral  techniques  for  bounding  sums: 


1 

V 


Econg(e)  =  ;^^E^' 

ee7  ^  i=l  l=i 


VV+l 


>2V 


(2VJ+1-2 

V 
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W  ith  some  straightforward  calculations,  the  reader  can  verify  that  the  term  in  parentheses  in  the  last  line  is 
greater  than  ^  for  all  j  and  k  such  that  k  >  j. 

The  calculation  for  the  upper  bound  is  similar: 


-  con 


j  k 


cC'- 


- 1  i 


2k^ 


3 

^Tk 


<  2k^ 


(  2  +  2y/r^ 

\ 


Thr  sum  in  parentheses  is  easily  shown  to  be  less  than  3  for  all  allowed  values  of  j  and  k  (i.e.,  k  >  j). 

The  conihination  of  these  bounds  implies  that,  for  every  path  in  the  embedding  that  starts  at  a  vertex 
Vj.  the  path  resistanre  is  0{k^).  Note  that  for  any  path  starting  at  one  of  the  leaves  this  value  is  increased 
by  at  most  1:  thus  any  path  in  the  embedding  has  resistance  proportional  to  k^.  Therefore  the  lower  bound 
on  A  is  0(p).  which  is  within  a  constant  of  the  upper  bound.  □ 

We  now  show  how  to  combine  the  star  embedding  lower  bounds  for  paths  with  masses  to  get  a  lower 
bound  on  the  uniform  priority  path  resistance  lower  bound  for  trees.  Consider  what  happens  when  we  pick 
any  point  along  critical  path  of  the  clique  embedding  and  set  it  to  zero.  The  zero  point  splits  the  path  in 
two.  producing  two  components  with  zero  Dirichlet  boundaries.  We  can  find  a  lower  bound  on  the  smallest 
eigenvalue  of  each  component. 

To  minimize  subscript  conflicts,  we  introduce  the  notation  r*  to  stand  for  We  have  three  r*’s 

to  consider:  one  from  the  clique  embedding  in  the  original  graph,  which  we  denote  as  r|;-,  and  the  two 
zero-boundary  bounds  with  respect  to  our  split  point.  Since  we  are  working  with  a  path,  we  can  think  of  the 
path  laid  out  with  vertices  in  increasing  order  from  left  to  right;  we  therefore  denote  the  two  path  resistances 
used  in  the  bounds  as  rj”  for  the  left  side  and  r*  for  the  right. 

Recall  that  the  lower  bound  for  the  clique  embedding  is  -T-;  we  can  relate  it  to  the  boundarv  case  as 
follows:  For  tree  T  with  a  clique  embedding  and  uniform  priorities,  let  P  be  the  critical  path  with  vertices 
indexed  from  1  to  A-  -f  1  and  let  be  the  path  resistance  of  P,  Let  1  <  .s  <  A’  +  1  be  the  index  of  a  split 
point.  Removing  separates  T  into  subgraphs;  let  T\  and  Tr  be  the  subtrees  containing  s  plus  the  vertices 
of  P  with  indices  less  than  ,s  and  greater  than  .s  respectively.  Let  r*  and  r*  be  the  maximum  path  resistances 
for  the  star  embeddings  of  T\  and  Tr  respectively  when  is  a  zero  boundary  point  and  uniform  priorities 
are  used  (if  either  of  these  trees  is  empty,  its  maximum  path  resistance  is  zero). 

Lemma  5.3.  For  the  siUiafion  dr.scrihcd  ahoDc, 


1  n 


Proof.  The  critical  path  has  A  edges  and  A  +  1  vertices,  with  edge  i  between  vertices  i  and  ?  +  1.  Since 
all  of  the  terms  involved  are  positive,  we  work  with  the  reciprocals  of  the  quantities  in  the  lemma  statement. 
The  reciprocal  of  the  left  term  can  be  written  as  follows: 

5-1  i  k  kn 

=  +  II  H 

i-  1  j-  }  i--s  j"  T  !■  1 
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We  can  write  the  reciprocal  of  the  right  expression  as 


A:+l 


m<j 


i=i  yj=i  j=i+i 

i=i  \j=i  ”  ^ 

s— 1  i  k  fc+1 

<EE™j'+E  E  ”‘r 


+E 

i=s  y 


m,' 


fc+1 


\  n 


nij 


i=l 


i=s  j=i+l 


To  see  the  final  inequality,  recall  that  n  —  proves  the  lemma.  D 

We  can  use  this  to  show  how  to  assemble  a  vector  that  gives  a  Rayleigh  quotient  for  the  original  Laplacian 
that  is  within  a  logdiam(r)  factor  of  A2.  This  is  shown  in  the  next  theorem: 

Theorem  5.4.  The  uniform  priority  path  resistance  method  produces  a  lower  bound  ^  <  X2  for  the 
Laplacian  of  any  unweighted  tree  T  that  is  off  by  a  factor  that  is  0(log  diam(T)). 

Proof  Let  L  be  the  Laplacian  of  T.  We  consider  splitting  the  critical  path  at  various  points.  Since  rf 
and  r*  depend  on  the  split  point  s,  we  make  that  clear  by  writing  them  as  rf  (5)  and  r*(s)  respectively  in  this 
proof.  Note  that  as  s  increases,  rf  (s)  increases  and  r*{s)  decreases.  We  also  have  that  rf(l)  =  r*(/c  +  l)  =  0. 
Let  3(s)  =  This  ratio  is  unbounded  when  s  =  1  and  decreases  as  s  increases,  reaching  0  when 

/  rj  (s) 

s  =  k^l.  We  consider  two  cases:  when  there  exists  an  5  such  that  2  >  /3{s)  >  and  when  there  is  no  such 
s. 

In  the  first  case,  Lemma  5.3  gives  us  the  following: 


(5.6) 


n 


rf(s)  +  r;(s)  rj 


<  A2. 


K 


Theorem  5.1  gives  a  way  to  construct  the  vectors  xi  and  Xp  so  that  the  respective  Rayleigh  quotients 
for  the  left  and  right  sides  are  O  ^max(iog(s-i),i)^  q  ^max(iog(fc^y-a),i)  j  respectively.  Recall  that  these 
vectors  are  positive.  The  left  Rayleigh  quotient  has  value 


RQ\ 


the  right  has  value 


RQr  = 


EiiVi  rriiX^ 


We  construct  a  vector  x  for  the  whole  tree  as  follows:  assign  the  values  from  xi  and  the  negatives  of 
the  values  of  Xr  to  the  corresponding  vertices  on  the  critical  path,  with  =  0.  Extend  the  vector  to  the 
rest  of  the  tree  by  assigning  vertices  off  the  critical  path  the  values  of  their  connection  points  as  described 
in  Section  4.3.  Using  the  Rayleigh  quotient  for  x  as  an  upper  bound  on  A2  requires  x’s  entries  to  sum  to  0. 
This  is  easily  accomplished  by  scaling  one  of  the  vectors  xi  or  Xp;  this  scaling  does  not  affect  the  respective 
Rayleigh  quotients  RQ\  and  RQr-  By  construction,  we  have  the  following  (remember  that  Xs  =0): 


12i=l  (^i  ~ 
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This  is  thr  sum  of  tlio  numerators  of  RQ\  and  RQ^  divided  by  the  sum  of  their  denominators.  It  is  well 
known  that,  for  positive  values  n.  b,  c,  and  d, 


Thus  we  have  that 


A2  <  max(7?(:?i,  RQr, ). 


By  (5.6).  (5.7).  and  the  eondition  on  /?(.s),  the  desired  result  holds  for  this  ca.se. 

Now  assume  that  there  is  no  .s  such  that  2  >  /7(.s)  >  i.  Then  there  is  some  s  such  that  P{s)  >  2  and 
3{s  1)  <  i.  By  Lemma  5..3,  both  r, •(...)  i).|r;(.,.!  i)  lower  bounds  on  and  hcncc  on 

A2.  The  second  of  these  plus  the  condition  on  ft{s  +  1)  imply  that  <  ^.  A  similar  argument  show.s 

that  is  also  a  lower  bound. 

We  eonstruet  a  vector  x  to  demonstrate  an  upper  bound.  The  construction  is  slightly  different  this  time: 
we  construct  X]  with  respect  to  the  boundary  being  at  vertex  .s  +  1,  and  x^  with  respect  to  the  boundary 
being  at  vertex  .s,  Corresponding  values  are  mapped  into  x  as  in  the  previous  ease.  Once  again,  we  ean  scale 
one  side  as  necessary  to  instire  that  the  values  in  x  sum  to  zero  without  affecting  the  left  or  right  Rayleigh 
quotirnts. 

Note  that  if  it  wore  the  case  that 


x^Lx  _  niim(i?.Qi)  -f  mim{RQr) 
x'^  X  dcnom(7?Qi)  +  dcnom(7?.Qr)  ’ 

thr  argument  as  in  the  previous  case  would  show  that  max(7?Qi, /t!Qr)  would  be  an  upper  bound  on  A2. 
This,  combined  with  the  lower  bounds  for  this  case  demonstrated  above,  would  be  sufficient  to  prove  that 
the  theorem  statement  holds  for  this  case.  However,  there  is  no  zero  point  in  this  case,  and  the  equality  does 
not  hold.  To  solve  this  problem,  we  will  show  that  in  fact  the  Rayleigh  quotient  is  no  more  than  a  factor  of 
two  larger  than  the  right  hand  side  of  the  inequality. 

To  see  this,  note  that  the  denominator  of  the  actual  Rayleigh  quotient  is  the  same  as  in  the  preceding 
expression.  The  numerator  has  the  following  changes:  we  lose  two  edges  from  the  ends  of  X]  and  Xr  to 
presumed  zero  points;  this  decreases  the  numerator  by  ^  We  replace  these  edges  by  an  edge 

between  and  ■.  1  that  contributes  {x^  —  i  1)^  <  2.t^  -f  2.r^.;  j.  Hence 


x'^  Lx  ^  2  (num(7?Qi)  4- num(7?.Qr)) 
x'^  X  denom(/?.Qi)  +  denom(7?,Qr) 


<  2  max(7?(5i,  /?.Qr). 


and  there  are  upper  and  lower  bounds  within  a  log  factor  of  the  diameter. 

This  proves  that  there  is  a  constant  such  that  the  theorem  holds  for  the  second  case,  which  implies  that 
the  theorem  holds  for  both  cases  using  the  larger  of  the  two  constants  from  the  cases.  □ 


6.  Open  Questions.  A  number  of  intere.sting  questions  remain  open.  How  good  in  general  is  the 
bound  based  on  uniform  priorities?  Is  there  an  easy  way  to  set  path  priorities  that  works  well  in  general? 
There  is  also  the  question  of  whether  there  is  an  easy  way  to  calculate  priorities  for  trees  in  both  the  Dirichlct 
case  and  the  case  without  boundary  conditions  that  improves  the  lower  bounds. 
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